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Abstract. We selected an unbiased, flux-limited and almost volume-complete sample of 13 distant. X-ray luminous (DXL, 
, z ~ 0.3) clusters and one supplementary cluster at z = 0.2578 from the REFLEX Survey (the REFLEX-DXL sample). We 

^ ■ performed a detailed study to explore their X-ray properties using XMM-Newton observations. Based on the precise radial 

i/^ ' distributions of the gas density and temperature, we obtained robust cluster masses and gas mass fractions. The average gas 

. mass fraction of the REFLEX-DXL sample at rsoo, 0.116 ± 0.007, agrees with the previous cluster studies and the WMAP 

' baryon fraction measurement. The scaled profiles of the surface brightness, temperature, entropy, gas mass and total mass are 

characterized by a self-similar behaviour at radii above 0.2-0.3 r^oo- The REFLEX-DXL sample confirms the previous studies 
■ of the normalization of the scaling relations (L-T, L-M, M-T and M„as~r) when the redshift evolution of the scaling relations 

is accounted for. We investigated the scatter of the scaling relations of the REFLEX-DXL sample. This gives the correlative 
scatter of (0.20,0.10) for variable of (M,T) of the Msm~T relation, for example. 

Oh' Key words. Cosmology: observations - Galaxies: clusters: general - X-rays: galaxies: clusters - (Cosmology:) dark matter 

o 
in 

c/3 1 . Introduction & Bryan 2001; Voit et al. 2002; Zhang & Wu 2003; Ponman et 

^ al. 2003). The X-ray properties of the most massive clusters are 

_ The number density of galaxy clusters probes the cosmic evo- ^^us well described in hierarchical modeling. Therefore, mas- 

"O ■ °f large-scale structure (LSS) and thus provides an ef- ^-^^ g^j^^y ^j^^^^^^ especially important in tracing LSS 

^ • fective test of cosmological models. It is sensitive to the matter evolution. The most massive clusters provide the cleanest re- 

^ density, Q^, and the amphtude of the cosmic power spectrum ^^j^^ compailng theory with observations. 
■ ■ ■ ■ on cluster scales, erg (e.g. Schuecker et al. 2003). Its evolution 

is sensitive to the dark energy, the density of which is charac- Excluding cooling cores (Fabian & Nulsen 1977), a self- 

terizedby the parameter and the equation of state parameter similar behaviour of the disti'ibutions of the ICM proper- 

w(z) (e.g. Vikhlinin et al. 2002; Allen et al. 2004; Chen & Ratra ties such as the temperature, density and entropy of mas- 

2004). sive clusters (> 4 keV) is indicated in ROSAT, ASCA, 

The most massive clusters show the strongest and clean- Chandra and XMM-Newton observations (e.g. Mai'kevitch 

est effects in the cosmological evolution. The structure of the 1998; Markevitch et al. 1998; Vikhlinin et al. 1999, 2005, 

X-ray emitting intra-cluster plasma in massive galaxy clusters 2006; Arnaud et al. 2002; Reiprich & Bohringer 2002; Zhang 

is predominantly determined by gravitational effects and shock et al. 2004a, 2004b, 2005b; Ota & Mitsuda 2005; Pratt & 

heating and is less affected by non-gravitational processes than Arnaud 2005; Pointecouteau et al. 2005) and simulations (e.g. 

low mass clusters. Only with decreasing cluster mass and intra- Borgani 2004; Borgani et al. 2004; Kay 2004; Kay et al. 2004). 

cluster medium (ICM) temperature, non-gravitational effects As a consequence of the similarity, massive galaxy clusters 

play an important role before and after the shock heating (Voit show tight scaling relations such as the luminosity-temperature 

(L-T, e.g. Isobe et al. 1990; Markevitch 1998; Arnaud & 

Send offprint requests to: Y.-Y. Zhang Evrard 1999; Ikebe et al. 2002; Reiprich & Bohringer 2002), 

* This work is based on observations made with the XMM-Newton, luminosity-mass (L-M, e.g. Reiprich & Bohringer 2002; 

an ESA science mission with instruments and contributions directly Popesso et al. 2005), mass-temperature (M-T, e.g. Nevalainen 

funded by ES A member states and the USA (NASA). et al. 2000; Finoguenov et al. 2001b; Neumann & Arnaud 
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2001; Xu et al. 2001; Horner 2001; Reiprich & Bohringer 
2002; Sanderson et al. 2003; Pierpaoli et al. 2001, 2003), 
and luminosity-metallicity (L-Z, e.g. Garnett 2002) relations. 
Additionally, a reliable estimate of the intrinsic scatter of the 
scaling relations is the key to a correct modeling to constrain 
cosmological parameters in the cosmological applications us- 
ing galaxy clusters. Therefore, understanding the intrinsic scat- 
ter of the scaling relations is as important as studying the scal- 
ing relations themselves. For example, on-going cluster merg- 
ers partially account for the scatter in the scaling relations since 
on-going mergers may lead to a temporary increase not only in 
the temperature and X-ray luminosity (Randall et al. 2002), but 
also in the core radius and slope of the surface brightness pro- 
file. Precise measurements of the ICM structure are required 
to allow accurate cluster mass and gas mass fraction determi- 
nations and thus to investigate the X-ray scaling relations for 
comparison of theory with simulations and observations. 

The ROSAT-ESO Flux-Limited X-ray (REFLEX, 
Bohringer et al. 2001, 2004) galaxy cluster survey pro- 
vides the largest homogeneously selected catalog of X-ray 
clusters of galaxies so far. It provides the basis to construct 
an unbiased sub-sample of clusters with specific selection 
criteria. We exploit it to compose a sample of distant. X-ray 
luminous (DXL) galaxy clusters in the redshift range, z - 0.27 
to 0.31, with Lx > 10'^^ erg s"' for the 0.1-2.4 keV band and 
one supplementary cluster at z = 0.2578 (the REFLEX-DXL 
samplejH The volume completeness correction can be done 
using the well known selection function of the REFLEX 
survey (Bohringer et al. 2004). 

Prime goals for the study of the REFLEX-DXL sample 
are, (1) to obtain reliable ICM properties such as temperature 
structure (Zhang et al. 2004a, Paper 11), (2) to determine accu- 
rate cluster masses and gas mass fractions using precise ICM 
property measurements, (3) to measure the normalization of 
the scaling relations with an improved accuracy and to dis- 
cuss their intrinsic scatter in detail, and (4) to test the evo- 
lution of both the scaling relations and the temperature func- 
tion (e.g. Henry 2004) by comparing the REFLEX-DXL sam- 
ple (at z ~ 0.3) to nearby cluster samples. We present the re- 
sults from high quality XMM-Newton data in this work. The 
data reduction is described in Sect. |2] We derive the X-ray 
properties of the ICM and determine the total masses and gas 
mass fractions based on precise gas density and temperature 
radial profiles in Sect. |3]l. We investigate the self-similarity of 
the REFLEX-DXL clusters (Sect. |4]l and discuss the peculiari- 
ties of the individual clusters accounting for the scatter around 
the self-similarity (Sect. [Sj. In Sect. |6] we draw conclusions. 
Unless explicitly stated otherwise, we adopt a flat ACDM cos- 
mology with the density parameter - 0.3 and the Hubble 
constant Hq - 70 km s"' Mpc"'. We adopt the solar abun- 
dance values of Anders & Grevesse (1989). Confidence in- 



' An Einstein-de Sitter cosmological model with fi,,, = 1.0, Q.\ = 
0.0 and the Hubble constant Hq = 50 km s"' Mpc"' was used for 
the Lx threshold in the sample construction. This luminosity threshold 
corresponds to Lx > 5.9 x lO"*"* erg s"' for a flat A cold dark matter 
(ACDM) cosmology with the density parameter = 0.3 and the 
Hubble constant Hg = 70 km s"' Mpc"' 



tervals coiTespond to the 68% confidence level. We apply the 
Orthogonal Distance Regression method (ODR; e.g. Feigelson 
& Babu 1992; Akritas & Bershady 1996) and take into account 
measurement errors on both variables for the parameter fitting 
of the scaling relations. We use Monte Carlo simulations for 
the uncertainty propagation on all quantities of interest. 

2. Data reduction 

2.1. Data preparation 

All 14 REFLEX-DXL clusters were observed by XMM- 
Newton in AO-1. Some properties of these observations and 
an overview of the sample are described in Bohringer et al. 
(2006, Paper 1). The observations of 5 clusters in AO-1 were 
heavily contaminated by flares. Zhang et al. (2004a) investi- 
gated the temperature structure of 9 remaining REFLEX-DXL 
clusters observed in AO-1. Finoguenov et al. (2005) applied a 
2-dimensional approach to study the structure such as projected 
density, temperature, pressure and entropy maps for those 9 
REFLEX-DXL clusters. All 5 clusters flared in AO-1 were re- 
observed in AO-3 in which four have sufficient quality for a 
detailed study. All the clusters of the REFLEX-DXL sample 
were uniformly analyzed in this work. 

All observations were performed with thin filter for three 
detectors. The MOS data were taken in Full Frame (FF) mode. 
The pn data were taken in Extended Full Frame (EFF) mode in 
AO-1 and FF mode in AO-3, respectively. For pn, the fractions 
of the out-of-time (OOT) effect are 2.32% and 6.30% for the 
EFF mode and FF mode, respectively. An OOT event file is 
created and used to statistically remove the OOT effect. 

Good calibration is required for correct data reduction. The 
difference of the spectral measurements using low energy cut- 
off values of 1 keV and 0.4 keV was significant at the 1-cr 
confidence level for most data sets in the previous analysis us- 
ing the calibration in XMMSAS v5.4.0 (Zhang et al. 2004a). 
Now we use the XMMSAS v6.5.0 software for the data re- 
duction. The new calibration in XMMSAS v6.5.0 provides a 
better agreement between data and model for EPIC, which are 
described in XMM-SOC-CAL-TN-0018 (Kirsch 2005). The 
CTI correction, point-spread function (PSF) core, astrometry, 
gain and energy redistribution have been improved. For ex- 
ample, in XMMSAS v6.0 the photon energy redistribution 
was strongly increased for energies below 0.53 keV in order 
to reduce excesses seen at low energies. Subsequent analysis 
of EPIC-pn spectra revealed large residuals around 0.43 keV 
which clearly showed that the redistribution was wrongly mod- 
eled in XMMSAS v6.0. The values in the redistribution ma- 
trix at low energies using the ground calibration were replaced 
based on the analysis of a set of spectra in the 0.1-2 keV 
band. For XMMSAS v6.5.0 the redistribution is performed 
energy-dependently to flattening the residuals. The spectra be- 
low 0.5 keV are greatly improved (Kirsch 2005). The spectral 
measurements using low energy cut-off values of 1 keV and 
0.4 keV are now consistent for most REFLEX-DXL clusters 
except for RXCJ0658-5556 (~ 15%). With the new calibration 
in XMMSAS v6.5.0, we therefore use the 0.4-10 keV band for 
the spectral analysis. 
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Above 10 keV (12 keV), there is little X-ray emission from 
clusters for MOS (pn) due to the low telescope efficiency at 
these energies. The particle background therefore dominates. 
The light curve in the range 10-12 keV (12-14 keV) for MOS 
(pn), binned in 100 s intervals, is used to monitor the particle 
background and to excise periods of high particle flux. Since 
episodes of "soft proton flares" (De Luca & Molendi 2004) 
were detected in the soft band, the light curve of the 0.3-10 keV 
band, binned in 10 s intervals, is used to monitor and to ex- 
cise the soft proton flares. A 10 s interval bin size is chosen 
for the soft band to provide a similar good photon statistic as 
for the hard band. The average and variance of the count rate 
(ctr) have been interactively determined for each light curve 
from the ctr histogram. Good time intervals (GTIs) are those 
intervals with ctrs below the threshold, which is defined as 
3-cr above the average. The GTIs of both the hard band and 
the soft band are used to screen the data. The background ob- 
servations are screened by the GTIs of the background data, 
which are produced using exactly the same thresholds as for 
the corresponding cluster target field. AU the observations of 
RXCJ201 1.3-5725 were almost completely contaminated by 
flares. In total, 183 pn counts (0.4-10 keV) are available for the 
imaging spectral analysis of this cluster RXCJ201 1.3-5725 is 
therefore excluded in the detailed study. Settings of FLAG - 
and PATTERN < 13 (PATTERN < 5) for MOS (pn) are used 
in the screening process. 



2.2. Source detection 



Over half of the clusters show clear substructures or/and elon- 
gation. An "edetect_chain" command has been used to detect 
point-like sources. Point sources except for those suspicious 
candidates described below are subtracted before the further 
data reduction. 

In the cluster center, point-like sources are hard to iden- 
tify because their emission is blended with the strongly peaked 
cluster emission, in particular in cooling core clusters (CCCs). 
At the angular resolution of the observations, it is difficult to 
distinguish between a steep cooling core cusp and a central 
point source (e.g. an AGN in a central dominant galaxy). In 
principle, cluster emission is thermal and point source emis- 
sion is non-thermal. The spectral shape can be used to identify 
point sources because a power law shape usually indicates the 
non-thermal emission contributed by point sources. We studied 
the spectra of the regions where there is a possible presence of 
point source in the cluster centers. We found that both "mekal" 
and "powerlaw" models provide acceptable fits. The photon 
statistic does not allow for a clear discrimination of the spectral 
shape. Therefore, a contribution of point sources to the cooling 
core cusp can not be ruled out for these clusters. The properties 
of these central regions are listed in Table [T] However, those 
suspicious candidates are not subtracted before the further data 
reduction. 



2.3. Background subtraction 

The background consists of several components exhibiting dif- 
ferent spectral and temporal characteristics (e.g. De Luca & 
Molendi 2001; Lumb et al. 2002; Read & Ponman 2003). The 
background components can approximately be divided into two 
groups (e.g. Zhang et al. 2004a). Group I contains the back- 
ground components showing significant vignetting (hereafter 
BVIGs), e.g. the cosmic X-ray background (CXB). Group II 
contains the components showing little or no vignetting (here- 
after BNVIGs), e.g. particle-induced background. 

It is safe to use a local background only if the vignetting 
effect is similar for the source and background regions. It 
becomes less reliable when large vignetting corrections are 
needed. The vignetting effect is energy-dependent and becomes 
significant above 5 keV. Additionally, the background com- 
ponent caused by the instrumental lines is different from po- 
sition to position. This cannot be accounted for by a local 
background. The REFLEX-DXL cluster emission covers most 
of the field of view (FOV), r ~ 8'. A local background is 
only available in the outskirts (e.g. 9.2' < r < 11.5') but re- 
quires substantial vignetting corrections. When the local back- 
ground is corrected for vignetting effects, the particle compo- 
nent is over corrected especially at high energies (artificially 
increased). This leads to an artificial softening of the derived 
cluster spectrum. The temperatures of the REFLEX-DXL clus- 
ters are above 5 keV, and thus tend to show a soft excess when 
the local background is used. Therefore, a local background is 
not the best choice for the analysis of the REFLEX-DXL sam- 
ple. 

Suitable background observations for such hot clusters are 
XMM-Newton observations of almost blank fields using the 
same instrumental set-up (e.g. a particular filter). We thus 
choose the blank sky accumulations in the Chandra Deep Field 
South (CDFS) as background. The CDFS observations used 
the same filter (thin filter for all detectors) and mode (FF/EFF 
mode for MOS/pn) as the REFLEX-DXL observations. In 
Zhang et al. (2005a), we investigated the blank sky accumu- 
lations in Lumb et al. (2002) as an alternative background and 
compared it to the CDFS observations. We found that the mea- 
surements are consistent within l-cr in both the spectral and the 
surface brightness analyses. 

Suitable background observations guarantee similar BVIGs 
as for the targets in the same detector coordinates. One can sub- 
tract such a background extracted in the same detector coordi- 
nates as for the targets. However, the difference between the tar- 
get and background should be taken into account as a residual 
background in the background subtraction. This can be done 
for the REFLEX-DXL clusters because the cluster X-ray emis- 
sion does not cover the whole FOV, but the r < 9' region. The 
outskirts in the FOV (e.g. 9.2' < r < 11.5') were used to ana- 
lyze the residual background. Strictly, the residual background, 
obtained after subtracting the blank field data, consists of both 
the residual CXB (because it varies on the sky) and the particle 
induced background (e.g residual background induced by soft 
protons). 
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2.4. Substructure 

As an unbiased sample with respect to cluster morphology, 
the REFLEX-DXL clusters are characterized in terms of their 
dynamical state (Table |2] also see Zhang et al. 2006) using 
the classification in Jones & Forman (1992). The substructure 
consideration is based on the cluster morphology classifica- 
tion. The detection of the cluster morphology and substructures 
strongly depends on the cluster brightness, spatial resolution, 
and exposure time (Jones & Forman 1992). We define the sub- 
structure regions using the 0.5-2 keV image as follows. 

As a start, we cut out the regions obviously contaminated 
by substructures. The remaining region can be fitted by a sym- 
metric surface brightness model convolved with the PSF plus 
the scaled CDFS background. We simulated a 0.5-2 keV image 
by convolving the symmetric surface brightness model with the 
PSF plus the scaled CDFS background. We subtracted the sim- 
ulated image from the observed 0.5-2 keV image and obtained 
a residual image. The summed counts of the substructure re- 
gion in the residual image characterize the substructure signal 
(S). The summed counts of the same region in the simulated 
image are used as the background (B). We scaled the residual 
image using the simulated image and defined a series of con- 
tour regions. The signal-to-noise-ratio for each contour region 
was calculated hy S /N - S / + B^. The region coiTespond- 
ing to the S /N >= 1 threshold is used to subtract the substruc- 
ture. After point source subtraction and substructure excision, 
hydrostatic equilibrium should usually be valid in the remain- 
ing sector (more details see Sect.|5]|. 

2.5. Spectral analysis 

We use the XSPEC vl 1.3.1 software for the spectral analysis. 
A double background subtraction procedure can be applied in 
two ways in the spectral analysis. One approach was firstly de- 
scribed in Zhang et al. (2004a, 2005a, hereafter DBS I). The 
other method was well illustrated in Pratt & Arnaud (2002, 
hereafter DBS II). We applied both methods to perform the 
data reduction and obtained similar results using these two ap- 
proaches. 

For a given region of interest covering cluster emission, the 
spectrum is extracted from the background observations in the 
same detector coordinates as for the target. Both the response 
matrix file (rmf) and auxiliary response file (arf) are used to re- 
cover the correct spectral shape and normalization of the cluster 
emission component. The following is usually taken into ac- 
count for the rmf and arf, (i) a pure redistribution matrix giving 
the probability that a photon of energy E, once detected, will be 
measured in data channel PI, (ii) the quantum efficiency (with- 
out any filter, which, in XMM-Newton calibration, is called 
closed filter position) of the CCD detector, (iii) filter transmis- 
sion, (iv) geometric factors such as bad pixel corrections and 
gap coiTections (e.g. around 4% for MOS), and (v) telescope 
effective area as a function of photon energy. The vignetting 
correction to eff'ective area for off-axis observations can be ac- 
counted for either in the arf or in the event lists by a weight 
column created by "evigweight". 



In DBS I, the rmf corresponds to (i) and (ii), and the arf 
corresponds to (iii), (iv), (v) as well as the vignetting correc- 
tion. A source spectrum is extracted from the outer region of 
the target observations. A background spectrum is extracted 
from the background observations in the same detector co- 
ordinates. Using only the rmf, the residual background spec- 
trum is found after subtracting the background spectrum from 
the source spectrum in the outer region. We assume as an ex- 
treme case that the residual background mainly consists of the 
BNVIGs, e.g. soft protons. The vignetting effect can thus be 
ignored for the residual spectrum. It is modeled by a "power- 
law/b" model ("powerlaw/b" ,a power law background model 
which is convolved with the instrumental redistribution ma- 
trix but not with the effective area) in XSPEC limited to 0.4- 
15 keV. This model is used to account for the residual back- 
ground in the spectral analysis over the whole energy range 
by applying a combined model of "wabs*mekal-i-powerlaw/b" 
("mekal", an emission model for hot diffiise gas, cf. Mewe et al. 
1985, 1986; Arnaud & Rothenflug 1985; Arnaud & Raymond 
1992; Kaastra 1992; Liedahl et al. 1995; "wabs", a model con- 
sidering the Galactic absorption) in XSPEC in the fitting pro- 
cedure. The coiTect shape of the background component is re- 
covered in the fit when the normalization of the "powerlaw/b" 
model is scaled to the area of the given region. The uncertain- 
ties of the fitting parameters of the "powerlaw/b" model are not 
introduced in the spectral analysis. This might result in an un- 
derestimate of the temperature uncertainty in the spectral anal- 
ysis using DBS I. However, the advantage of DBS I is that the 
shape of the residual spectrum is conserved during the proce- 
dure. 

In DBS II, all spectra are extracted considering the vi- 
gnetting correction by a weighted column in the event list pro- 
duced by "evigweight". The on-axis rmf and arf are co-created 
to account for (i) to (v). The target spectrum is extracted from 
the region of interest. The first-order background spectrum is 
given by a spectrum from the background observations which 
is extracted in the same detector coordinates as for the source 
spectrum, and which is scaled using the ctr ratio of the tar- 
get and background limited to 10-12 keV (12-14 keV) for 
MOS (pn). The second-order background spectrum, the resid- 
ual background spectrum, is prepared as follows. A source 
spectrum is extracted from the outer region of the target ob- 
servations and its background spectrum from the same outer 
region of the background observations. The background spec- 
trum is scaled using the ctr ratio of the target and background 
Hmited to 10-12 keV (12-14 keV) for MOS (pn), and is sub- 
tracted from the source spectrum to obtain a residual back- 
ground spectrum. This residual background spectrum is nor- 
malized to the area of the region of interest as a second-order 
background spectrum. Both the first-order and second-order 
background spectra are subtracted from the target spectrum. A 
combined model of "wabs*mekar' is then used with the on-axis 
arf and rmf in XSPEC for the fitting (Fig. [B- 

Basically, the vignetting correction is applied on the resid- 
ual background in DBS II and not applied in DBS I. The for- 
mer method is thus in principle better when the residual is 
dominated by the BVIGs (e.g. the CXB residual component), 
while the latter is better when the residual is dominated by 
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the BNVIGs. The two approaches provide consistent results at 
the 1-cr confidence level which indirectly tested the two ap- 
proaches. A small discrepancy was found mostly in the last 
bin but showed no trend of being higher for one particular 
approach. There is no preferred indication of one of the two 
approaches. Additionally, the DBS II approach is simple in, 
(i) modeling the residual background, by applying the residual 
spectrum, instead of looking for an acceptable (x^ < 2) power 
law model fit; (ii) allowing negative residual components, and 
(iii) generating one on-axis arf working for all the spectra ex- 
tracted from different annuli. Therefore, we adopt the measure- 
ments using DBS II to illustrate the further computations. 

2.6. Image analysis 

The 0.5-2 keV band is used to derive the surface brightness 
profiles (also see Zhang et al. 2005a). This ensures an almost 
temperature-independent X-ray emission coefficient over the 
expected temperature range. The vignetting correction to effec- 
tive area is accounted for using a weight column in the event 
lists created by "evigweight". Geometric factors such as bad 
pixel corrections are accounted for in the exposure maps. The 
width of the radial bins is 2". An azimuthally averaged sur- 
face brightness profile of CDFS is derived in the same detector 
coordinates as for the targets. The ctr ratios of the targets and 
CDFS in the 10-12 keV band and 12-14 keV band for MOS 
and pn, respectively, are used to scale the CDFS surface bright- 
ness. The residual background in the 0.5-2 keV band is intro- 
duced by using the determined residual spectrum in the spectral 
analysis. The background-subtracted and vignetting-corrected 
surface profiles for three detectors are added into a single pro- 
file, and re-binned to reach a significant level of at least 3-cr in 
each annulus. We take into account a 10% uncertainty of the 
scaled CDFS background and residual background. 

2. 7. PSF and de-projection 

Using the XMM-Newton PSF calibrations by Ghizzardi (2001) 
we estimated the redistribution fraction of the flux. We found 
20% for bins with width about 0.5' and less than 10% for bins 
with width greater than 1' neglecting energy dependent effects. 
The PSF blurring can not be completely considered for the 
spectral analysis as it is done for the image analysis because of 
the limited photon statistic. For such distant clusters, the PSF 
effect is only important within O.Srjoo and introduces an added 
uncertainty to the final results of the temperature profiles. This 
has to be investigated using deeper exposures with better pho- 
ton statistic. 

The projected temperature is the observed temperature 
from a particular annulus, containing in projection the fore- 
ground and background temperature measurements. Under the 
assumption of spherical symmetry, the gas temperature in each 
spherical shell is derived by de-projecting the projected spec- 
tra. In this procedure, the inner shells contribute nothing to the 
outer annuli. The projected spectrum in the outermost annu- 
lus is thus equal to the spectrum in the outermost shell. The 
projected spectrum in the neighboring inner annulus has con- 



tributions from all the spectra in the shell at the radius of this 
annulus and in the outer shells as shown in Suto et al. (1998). In 
XSPEC, a combined fit to the projected spectra measured in all 
annuli simultaneously gives the contribution of each spherical 
shell to each annulus and provides the de-projected measure- 
ments of the temperature and metallicity. 

In the imaging analysis, we correct the PSF effect by fit- 
ting the observational surface brightness profile with a surface 
brightness model convolved with the empirical PSF matrices 
(Ghizzardi 2001). The surface brightness model is calculated 
by the projection of the radial electron density profile. 

3. X-ray properties 

The primary parameters of the REFLEX-DXL clusters are 
given in Table [3l 

3.1. Density contrast 

The mean cluster density contrast. A, is the average density 
with respect to the critical density, pdz) = pcoE^iz), where 
E\z) = Qm(l -H z)^ -H Qa -H (1 - Q„ - ^aK^ + z)\ r„ is the 
radius within which the density contrast is A. Ma is the total 
mass within ta. For A - 500, rsoo is the radius within which 
the density contrast is 500 and M500 is the total mass within 

'"500 ■ 

3.2. Metallicity and temperature 

We found that the X-ray determined redshifts agree with the 
optically measured redshifts (also in Zhang et al. 2004a). We 
therefore fixed the redshift to the optical redshift in the further 
analysis. The spectrum is fitted by a combined "wabs*mekal" 
model with fixed Galactic absorption (Dickey & Lockman 
1990) and redshift (Bohringer et al. 2004a). 

Temperature profiles can provide a useful means to study 
the thermodynamic al history of galaxy clusters. XMM-Newton 
(and Chandra), in contrast to earlier telescopes, provides a less 
energy-dependent, and smaller, PSF. It is thus more reliable 
to study cluster temperature profiles with XMM-Newton. A 
systematic spectral analysis was performed in annuli. We re- 
binned the data to contain around 500-550 net counts per an- 
nulus for MOSl in the 2-7 keV band. In XSPEC, a combined 
fit of "wabs*mekal*projcf 'to the projected spectra measured in 
all annuli simultaneously gives the contribution of each spher- 
ical shell to each annulus and provides the de-projected mea- 
surements of the temperature and also metallicity. In Fig.|2]we 
show the de-projected temperature profiles of the REFLEX- 
DXL clusters. The temperature profiles are approximated by 
the parameterization T(r) - r3exp[-(r - T\f-IT2\ + 7^6(1 -H 
p-/T^y'^^ + T-i crossing all the data points. The temperature 
measure uncertainties are approximated by T(r) x (Tg + Tgr) 
and are propagated for individual data points by Monte Carlo 
simulations in the further analysis. 

Vikhlinin et al. (2005) used the weighting of the 3-D tem- 
perature profile and averaged the temperature profile in a cer- 
tain radial range. For the REFLEX-DXL sample, the radial 
temperature profiles are almost self-similar above 0.1-0.2r5oo- 
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The REFLEX-DXL data set shows significantly low S /N above 
O.Srsoo- Therefore, the volume averaged radial temperature pro- 
file of 0.1-0.5 r5()o is chosen as the global temperature T listed 
in Table [3] Similarly, the metallicity profile of 0.1-0.5 rsoo is 
used to obtain the global metallicity Z also listed in Table[3] We 
found an average of 0.24+0.03 Zq for the global metallicities of 
the REFLEX-DXL sample excluding RXCJ201 1.3-5725. This 
agrees with the value of 0.21 Zq in Allen & Fabian (1998). It 
is also consistent with the averaged metallicity Z - 0.21^q q"Zq 
for 18 distant (0.3 < z < 1.3) clusters in Tozzi et al. (2003).^ 

Based on the flared observations of RXCJ201 1.3-5725, we 
measured a global spectral temperature (~ 4 + 2 keV) with a 
fixed metallicity of 0.3 Z© within 0.5' radius. 

3.3. Surface brightness 

A y6-model (e.g. Cavaliere & Fusco-Femiano 1976; Jones & 
Forman 1984) is often used to describe electron density pro- 
files. To obtain an acceptable fit for all clusters in this sam- 
ple, we adopt a double-j6 model, rigir) = neOi(l + I r^iT^^^'^ + 
neoiil / r^2^^^^^^ ■ The X-ray surface brightness profile model 
is linked to the radial profile of the ICM electron number den- 
sity as an integral performed along the line-of-sight, 

SxiR)'^ nlM. (1) 
Jr 

We fit the observed surface brightness profile by this integral 
convolved with the PSF matrices (Fig. O and obtain the pa- 
rameters of the double-yS model of the electron density pro- 
file. The fit was performed within the truncation radius (rt, see 
Table [3l) corresponding to a 5 /A^ of 3 of the observational sur- 
face brightness profile. The truncation radii, rt, are about or 
above rsoo for the REFLEX-DXL clusters. The X-ray bolomet- 
ric luminosity (here we use the 0.01-100 keV band) is given 
by L^°' oc ^ K{r)n^{r)dV, practically an integral of the X-ray 
surface brightness to infinity (In practice, 20' is used for the in- 
tegral upper radial limit). We applied a power law to fit the sur- 
face brightness in the outskirts (> 0.4r5oo). More details about 
the slope of the the outskirts can be found in Sect. 15.11 

3.4. Cooling time 

The cooling time is derived by the total energy of the gas di- 
vided by the energy loss rate 

tc = — ^ (2) 
niA 

where A, rig, and T are the radiative cooling function, gas 
number density, electron number density and temperature, re- 
spectively. We compute the upper limit of the age of the clus- 
ter as an integral from the cluster redshift z up to z = 100. 
Cooling regions are those showing cooling time less than the 
upper limit of the cluster age. The boundary radius of such a 
region is called the cooling radius. The cooling radius is zero 
when the cooling time is larger than the cluster age. The cool- 
ing times and cooling radii are given in Table |4] 



3.5. Gas entropy 

The entropy is the key to an understanding of the thermody- 
namical history since it results from shock heating of the gas 
during cluster formation. The observed entropy is generally de- 
fined as S - Tn^^^^ for galaxy cluster studies (e.g. Ponman et 
al. 1999), and it scales with the cluster temperature. An excess 
above the scaling law indicates non-gravitational heating ef- 
fects by sources such as the feedback from super-novae (SN) 
and AGNs (e.g. Lloyd-Davies et al. 2000). Radiative cooling 
can also raise the entropy of the ICM (e.g Pearce et al. 2000) 
or produce a deficit below the scaling law. In this sample, the 
clusters appearing more relaxed are identical to those show- 
ing lower central entropy values. For the REFLEX-DXL clus- 
ters, the most inner temperature data points are measured at 
~ O.lrsoo. This leads to the fact that the entropies at O.lrsoo 
(5^ O.I 1500 ) ^'"^ resolved and can be used as the values of the 
central entropies as shown in the entropy radial distributions 
(Fig. [Toll. At 0.1r2oo (about 0.18r50()), the entropy increases as a 
function of radius and is thus significantly larger than the cen- 
tral entropy. 

3.6. Mass distribution 

We assume that the intra-cluster gas is in hydrostatic equilib- 
rium within the gravitational potential dominated by DM. The 
ICM can thus be used to trace the cluster potential. Under the 
assumption of spherical symmetry, the cluster mass is calcu- 
lated from the X-ray measured ICM density and temperature 
distributions, 

1 d[n,(r)T(r)] _ GM(< r) 

where yu - 0.62 is the mean molecular weight per hydro- 
gen atom. Following the Monte Carlo simulation method (e.g. 
Neumann & Bohringer 1995), we use a set of input parameters 
of the approximation functions, in which /3, WeOi, (/ = 1 , 2) 
are for the gas density radial profile ne('") and (/ = 1, ...,7) 
are for the temperature radial profile T{r), respectively, to com- 
pute the cluster mass. The individual data point uncertainties 
are propagated by Monte Carlo simulations. We used the mea- 
sured mass profile to estimate Msoq and r^oo. 

The NEW model (e.g. Navarro et al. 1997, 2004, NEW) 
cannot provide an acceptable fit for the observed mass pro- 
files. We therefore adopt the best fit of an extended-NFW 
model (e.g. Hernquist 1990; Zhao 1996; Moore et al. 1999), 
PDuir) = ps('"/'"s)""(l + r/r^)"'^, where ps and are the 
the characteristic density and scale of the halo, respectively. 
In Fig. |5] we show the measured mass profiles and their best 
extended-NFW model fits. We derive the concentration param- 
eters (cgoo) of the DM distributions of the REFLEX-DXL clus- 
ters (see Table [Sj using the extended-NFW model. The dif- 
ference in the concentration parameters using different mod- 
els can be large. For instance, an extended-NFW profile fit- 
ted to the numerical simulations of Moore et al (1999) gives 
a concentration parameter 50% higher than the concentration 
parameter given by the NFW model. In general, the theoretical 
and observed concentration-mass relations are compared when 
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the same models (e.g. NFW) are considered. However, the 
REFLEX-DXL data cannot be well fitted by the NFW model. 
We thus cannot directly compare our results to the published 
observations (e.g. Pointecouteau et al 2005) and simulations 
(e.g. Dolag et al. 2004) which are based on the NFW model. 
The best fit of the REFLEX-DXL sample gives the slope of 
1.5 ± 0.2 and the normalization of 10' '=^"^ (Fig.|6]). If we fixed 
the slope parameter to -0.102 for the concentration-mass rela- 
tion as found in Dolag et al. (2004), the best fit of the REFLEX- 
DXL sample gives the normalization of 8.7 + 7.8. 

3. 7. Gas mass fraction distribution 

The gas mass fraction is an important parameter for cluster 
physics, e.g. heating and cooling processes, and cosmological 
applications using galaxy clusters (e.g. Vikhlinin et al. 2002; 
Allen et al. 2004). The gas mass fraction distribution is defined 
to be /gas(< r) = Mgas(< r)/M(< r). 

As shown in Fig. |7] the gas mass fractions increase as a 
function of radius in the rsooo < r < riooo region. This indicates 
that the DM distribution is more concentrated in the center than 
the gas distribution. We derived an average gas mass fraction 
of 0.116 ± 0.007 at rsoo. This agrees with the gas mass frac- 
tions found for many massive clusters showing temperatures 
greater than 5 keV (e.g. Mohr et al. 1999). The gas mass frac- 
tions around r25oo show the smallest scatter, 0.100 + 0.007, and 
the values are similar to the measurements of Allen et al. (2002) 
based on Chandra observations of 7 clusters yielding /gas ~ 
0.105-0.138/!;o^^. At rim, the extrapolated gas mass fractions 
show consistency with the measurements of Sanderson et al. 
(2003) based on ASCA/GIS, ASCA/SIS and ROSAT/PSPC ob- 
servations of 66 clusters yielding /gas - 0.13 + O.OI/z^q , the 
measurements of Ettori et al. (2002a) based on BeppoSAX 
observations of 22 nearby clusters, and the gas mass fraction 
for A1413 (Pratt & Arnaud 2002) at z = 0.143 based on 
XMM-Newton observations yielding /gas ~ QAlh^^J^. As ex- 
pected, the gas mass fraction distributions of the REFLEX- 
DXL clusters are lower than the universal baryon fraction, 
/b = Ob/Om = 0.167 ± 0.014, based on the recent WMAP 
measurements, h = 0.71+°{;4^ D.^ = 0.022 ± 0.001 and 
D.^ = 0.132+J}J»^ (e.g. Spergel et al. 2003; Hansen et al. 
2004). This is because the baryons in galaxy clusters reside 
mostly in hot gas together with a fraction (15% /gas) in stars 
as implied from simulations (e.g. Eke et al. 1998; Kravtsov et 
al. 2005). In principle, Q^, can be determined from the baryon 
fraction, fi, - /gas + /gai, in which a contribution from stars 
in galaxies is given by /gai - 0.02 ± 0.0 1/;^,)' (White et al. 
1993). The gas mass fractions, ~ 0.116 ± 0.007, at rsoo of the 
REFLEX-DXL clusters support a low matter density Universe 
as also shown in recent studies (e.g. Allen et al. 2002; Ettori et 
al. 2003; Vikhhnin et al. 2003). 

4. Scaling relations 

Simulations (e.g. Navarro et al. 1997, 2004) suggest a self- 
similar structure for galaxy clusters in hierarchical structure 
formation scenarios. The scaled profiles of the X-ray proper- 
ties and their scatter can be used to quantify the structural vari- 



ations. This is a probe to test the regularity of galaxy clusters 
and to understand their formation and evolution. The accuracy 
of the determination of the scaling relations, limited by how 
precise the cluster mass and other global observable can be es- 
timated, is of prime importance for the cosmological applica- 
tions of clusters of galaxies. 

Because the observational truncation radii (rt) in the surface 
brightness profiles are about or slightly above rsoo, we use rjoo 
for radial scaling. 

The following redshift evolution corrections (e.g. Ettori et 
al. 2004) are usually used to check the dependence on the evo- 
lution of the cosmological parameters, 

Sx-E-\z)(AcJAc,or''<^f(T), 

S ■ E^l\z) {AcJKof^ ^ f(T), 

L-£-i(z)(Ac,JA,,o)-°-5cx/(r), 

M ■ E(z) (AcJAcfif-' 'X f(T), 

Mgas ■ Eiz) (A,,,/A,,o)°-5 cc f{T), 

where Ac,z = ISn^ + S2[Q.^_^ - 1] - 39[Qn,,z - 1]^ for a flat 
Universe and Q^.z is the cosmic density parameter at redshift 
z- 

4.1. Scaled temperature profiles 

Studies of the cluster temperature distributions (e.g. 
Markevitch et al. 1998; De Grandi & Molendi 2002) in- 
dicate a steep decline beyond an isothermal center We scaled 
the radial temperature profiles by the global temperature T 
and rsoo as shown in Fig. [8] An average temperature profile 
is derived by averaging the 1-cr boundary of the scaled radial 
temperature profiles of the REFLEX-DXL clusters. As shown 
in Fig. |8] we found a closely self-similar behaviour. Up to 
0.3r5()(), we observed an almost constant temperature distribu- 
tion with a temperature peak at around 0.2 rsoo. Four clusters 
(RXCJ0232.2-4420, RXCJ0307.0-2840, RXCJ0437.1-H0043 
and RXCJ0528. 9-3927) show relatively cool cluster cores. 
A temperature profile decreasing down to 80% of the peak 
value has been observed outside 0.3r5oo for the average of the 
REFLEX-DXL clusters. Three clusters (RXCJ0437. 1-1-0043, 
RXCJ0658.5+5556, and RXCJ2308.3+0211) show tempera- 
tures rising with relatively large error bars in their temperature 
profiles. The average temperature profile can decline down to 
50% of the peak value when these three clusters are excluded. 
This average temperature profile is consistent with the average 
profiles from ASCA in Markevitch et al. (1998), BeppoSAX in 
De Grandi & Molendi (2002) and Chandra (using an assumed 
uncertainty of 20% of the averaged temperature profile as an 
approximate illustration) in Vikhlinin et al. (2005) within the 
observational dispersion. A similarly universal temperature 
profile is indicated by simulations (e.g. Borgani et al. 2004; 
Borgani 2004). 

4.2. Scaled surface brightness profiles 

We scaled the surface brightness profiles using the self-similar 
scaling, Sx °^ T^'^, and the empirical scaling, Sx °^ T^^^^, 
respectively, as described in e.g. Arnaud et al. (2002). Both 
scaled profiles show a less scattered self-similar behaviour at 
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r > Q.lrsoo (see Fig. |9]l. The core radii populate a broad range 
of values, 0.03-0.2 rjoo. 



4.3. Scaled entropy profiles 

According to the standard self-similar model the entropy scales 
as 5 oc r (e.g. Pratt & Arnaud 2005). Ponman et al. (2003) 
suggested to scale the entropy as 5 oc t'O-Ss \Ye investi- 
gate the entropy-temperature relation (S-T) using 5o.ii-5o„ to 
represent the central entropies. Four REFLEX-DXL clusters 
(RXCJ0232.2-4420, RXCJ0307.0-2840, RXCJ0437. 1+0043 
and RXCJ0528. 9-3927) show significantly lower central 
entropies compared to the S-T scaling law. These clus- 
ters have relatively cool cluster cores as observed in their 
temperature profiles. Neglecting the resolution problem, we 
found that the radiative cooling effect introduces a signif- 
icant scatter in the S-T scaling relation in terms of a 
lower (~ 7%) normalization. Furthermore, 3 pronounced 
merger clusters (RXCJ0014.3-3022, RXCJ05 16.7-5430 and 
RXCJ2337. 6-1-0016) show relatively higher central entropies. 
This indicates that mergers may also introduce some scatter in 
the S-T relation but in terms of a slightly higher (~ 10%) nor- 
malization. The central entropy can thus be used not only as a 
mechanical educt of the non-gravitational process, but also as 
an indicator of the merger stage. Excluding the above 7 clusters 
discussed, we performed a best fit for the S-T relation of the 
remaining 6 clusters giving S oc 70.63±o.i5 -pj^jg ^^j. ^gj-ggg ^^j^jj 
the S-T relation of the Birmingham-CfA clusters (Ponman et 
al. 2003). Therefore, we scaled the radial entropy profile us- 
ing the empirically determined scaling (Ponman et al. 2003), 
S oc 70 65^ ^jjj ^^^^y shown in Fig. [10] the scaled entropy 
profiles of the REFLEX-DXL clusters agree with the scaled 
entropy profiles of the Birmingham-CfA clusters in Ponman et 
al. (2003) and the clusters in Pratt et al. (2006) in the same 
temperature range (6-20 keV) within the observational dis- 
persions. The least scatter (~ 30%) of the entropy profiles is 
around 0.2-0.3r5oo. The combined entropy profiles give the 
best fit, S{r) oc ^o.99±o.o6^ above 0.1r5oo. A similar power law 
as 5 oc r"-^^ was found by Ettori et al. (2002b) and S oc r" ''^ by 
Piffaretti et al. (2005). The slope of the entropy profiles is shal- 
lower at the 1-0" confidence level than the predicted slope from 
a spherical accretion shock model, S oc r' ' (e.g. Kay 2004). 



4.4. Scaled total mass and gas mass profiles 

The mass profiles were scaled with respect to M500 and rsoo, 
respectively (Fig.fTTTi. We found the least scatter at radii above 
0.2-0.3 rg(K). In the inner parts (< 0.2r5oo), the mass profiles 
vary significantly with the cluster central dynamics (Fig. fTTT l. 
For the clusters showing cooling cores, the mass distributions 
are relatively cuspidal, while for the merger clusters, less con- 
centrated mass distributions are observed. The scaled gas mass 
profiles appear more self-similar than the scaled mass profiles, 
especially at radii above 0.2-0.3 rgoo- 



4.5. ROSAT and XMM-Newton luminosities 

Substructure is often observed in galaxy clusters and the fre- 
quency of its occurrence has for example been estimated to be 
of the order of about 52 + 7% (Schuecker et al. 2001). The high 
resolution XMM-Newton data allow us to identify the substruc- 
tures and point-like sources better than what was possible with 
the earlier X-ray telescopes. Subtracting the substructure con- 
tribution also from the ROSAT measured cluster luminosity, 
we found a good agreement between the ROSAT and XMM- 
Newton luminosity for the REFLEX-DXL clusters (Fig. [T2T i. 
The XMM-Newton data provide a more reliable and complete 
detection of substructures and point-like sources. Therefore we 
make use of this capability to obtain a better approximation 
to spherical symmetry and dynamical equilibrium of the main, 
largely relaxed component by excluding the substructures and 
point-like sources. The properties of the main cluster compo- 
nent are more representative for investigating various scaling 
relations of regular galaxy clusters, such as the L-F relation. 

4.6. Scaling relations 

To use the temperature/mass function of this unbiased, flux- 
limited and almost volume-complete sample to constrain cos- 
mological parameters, it is important to calibrate the scaling 
relations between the X-ray luminosity, temperature and grav- 
itational mass. The scaling relations can generally be parame- 
terized by a power law. 

The scatter describes the dispersion between the observa- 
tional data points and the best fit. The scatter in the scaling 
relation is strongly dependent on the temperature measurement 
uncertainty. The massive clusters in a narrow temperature range 
provide an important means to constrain the normalization of 
the scaling relations such as the M-T relation. We collected 
recently published scaling relations and compared them to the 
results of this work in Table |6] and Figs. [T3l - [T9l The scatter in 
the scaling relations can partially be explained by variation of 
the cluster morphology. Comparing the REFLEX-DXL sample 
to the nearby samples, we found that the evolution of the scal- 
ing relations are accounted for by the redshift evolution given 
in Sect. H) We obtained an overall agreement with the recent 
studies of the scaling relations within the observational disper- 
sion. This fits into the general opinion that galaxy clusters are 
self-similar up to z ~ 1 (e.g. Arnaud 2005). 

To determine the normalization, the slope of the M^oq-T 
relation is fixed to 1.5 in the fitting procedure as also used 
in many published results (e.g. Evrard et al. 1996, Ettori et 
al. 2004). The M^oq-T relation (Fig. [T3]) agrees with those in 
Evrard et al. (1996), Ettori et al. (2004) and Vikhlinin et al 
(2006). As an unbiased sample, the M^qq-T of the REFLEX- 
DXL sample shows slightly higher normalization than the local 
relations found in Finoguenov et al. (2001b) and Popesso et al. 
(2005), and the local relation derived in Arnaud et al (2005) 
for relaxed clusters also based on the XMM-Newton tempera- 
ture profiles. Otherwise, there is no obvious additional evolu- 
tion of the normalization of the M-T relation after the redshift 
evolution correction. This was also found for the other distant 
clusters in Maughan et al. (2003) and Ettori et al. (2003). 
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As shown in Fig. [T] the gas approximately follows DM at 
radii above ^2500 for most of the REFLEX-DXL clusters. The 
gas mass can thus be used as a measure of the total mass. We 
fixed the slope parameter to 1 .8 as also used in Castillo-Morales 
& Schindler (2003) and Borgani et al. (2004) and fitted the 
normalization. The Mgas.soo-?^ relation for the REFLEX-DXL 
clusters (see Table |6] and Fig. [T4l i is in good agreement with 
the relation in Mohr et al. (1999) also using an X-ray flux lim- 
ited cluster sample. Our result also agrees with the recent result 
in Castillo-Morales & Schindler (2003). Recent simulations 
also indicate a strong M^^s-T scaling relation (e.g. Borgani 
et al. 2004). The normalization of the Mc,as,500-r relations is 
slightly higher for the observations than for the simulations 
(e.g. Borgani et al. 2004). 

The X-ray luminosity is a key parameter among the funda- 
mental cluster properties including also mass, temperature, and 
velocity dispersion. Excluding cooling cores (~ O.I/z^q Mpc), 
Markevitch (1998) reduced the scatter in the L^°'-r relation. 
As listed in Table |3] we list the bolometric X-ray luminosity of 
the REFLEX-DXL clusters including and excluding the cluster 
cores, r < O.lrgoo, as is used in many studies (e.g. Markevitch 
1998; Zhang 2001). The scatter of the L-T relation is reduced 
by 15% excluding cooling cores. About 8-33% of the lumi- 
nosity is contributed by the cluster cores (< O.lrgoo) for the 
REFLEX-DXL clusters. Therefore, the normalization of the L- 
T relation excluding cooling cores is also reduced by 10% for 
the REFLEX-DXL clusters. As listed in Table |6]and shown in 
Figs. [TsUTSl we fixed the slope parameters and fitted the nor- 
malizaton for the REFLEX-DXL sample after the redshift evo- 
lution correction. Within the observational scatter, the L^°'-r 
relation for the REFLEX-DXL sample (Fig. [TSI l agrees with 
the relation in Reiprich & Bohringer (2002) also as an unbi- 
ased sample, and the relations in Arnaud & Evrard (1999) and 
Markevitch (1998). An alternative redshift evolution in the L- 
T relation is described in Kotov & Vikhlinin (2005) yielding 
L cc 7"2 64^ +2)'**. The result here agrees with theirs when the 
alternative redshift evolution is adopted. 

In Fig.[T7]and Fig. [18] we show the L-M relations (Table|6]l 
of the REFLEX-DXL sample. The slope parameter is fixed to 
1.3, as derived in Popesso et al. (2005), in the fitting proce- 
dure. The best fit of the normaUzation of the L^^^^ '^^'^^-M re- 
lation for the REFLEX-DXL sample agrees with the best fits in 
Reiprich & Bohringer (2002) and Popesso et al. (2005). 

In Fig. [19] we show the gas mass fraction, /gas.soo, and gas 
mass fraction ratio, /gas,50o//eas,2500, as a function of cluster 
temperature. We found a weak evidence that the gas mass frac- 
tion increases with the cluster temperature, /gas,50o ^ T^O-^iO-S 
This agrees with the scaling, /gas.soo °^ 70-34±o.22^ found in Mohr 
et al (1999). The ratio of the gas mass fractions at larger and 
smaller radii (e.g. r^oo and rasoo) can be used to characterize 
the extent of gas relative to DM. When this ratio is greater 
than 1, the gas is more extended than DM (Reiprich 2001). For 
the REFLEX-DXL clusters, the gas mass ratios (Fig.[T9]l show 
that gas is more extended than DM in the cluster inner region 
but follows DM better in the outer region. Fig. [19] indicates a 
small trend that the gas is more extended than DM in massive 
clusters. This indicates that the ICM is less influenced by non- 
gravitational efl'ects and that the energy input is less important 



in the outer region for such massive clusters. This is also indi- 
cated by simulations (Rowley et al. 2004). 

4.7. Intrinsic scatter of the scaling relations 

The key to extracting cosmological parameters from the num- 
ber density of galaxy clusters is a correct understanding of the 
mass-observable scaling relations and their intrinsic scatter. 
The scatter of the mass-observable scaling relation describes 
how well the observable can be used as an estimator of the 
total mass. The coiTelative scatter (o^coi) includes the intrinsic 
scatter (o^int) and observational scatter (o^obs)- We investigated 
the logarithmic intrinsic scatter of the scaling relations, L-T, 
L-M, M-T, and M^^.-T, of the REFLEX-DXL sample. 

The observational scatter o^obs is the average of the esti- 
mated observational uncertainties CTobs- The correlative scatter 
(Tcor is the average of the deviation cTcor of the observational 
data points from the best fit of the scaling relation. Assuming 
that cTobii and 0-^01 are not correlated, we apply a Gaussian sta- 
tistical addition of the two eff'ects to compute the logarithmic 

intrinsic scatter (Tint, the average of 0-^1 - -^oIot^"^-^ . The 
reason that we used the average instead of the mean to cal- 
culate the scatter lies as following. We investigated the actual 
distributions of the observational scatter, coiTelative scatter and 
logarithmic intrinsic scatter. We found that in most cases the 
distribution deviates from a Gaussian. Fig. [20] shows for ex- 
ample the distribution of the logarithmic intrinsic scatter of 
A^gas.500 for the Mgas.soo-T" relation. The merger clusters (e.g. 
RXCJ05 16.7-5430 and RXCJ0658.5-5556) often faU into the 
right tail of the histogram. The asymmetry of the histogram 
is most probably due to the variety of cluster morphologies, 
e.g. mergers, which can produce significant deviation from the 
mean. It could also be due to the fact that the REFLEX-DXL 
sample does not contain enough members to give a pronounced 
Gaussian statistics. We thus use the average of the histogram to 
derive the logarithmic intrinsic scatter for the whole sample, 
which is listed in Table|2]for the scaling relations. 

The mass-observable scaling relation which shows the 
least scatter provides the best prediction for the total mass or 
gas mass. We find here that the temperature is the best esti- 
mator However, we have to recall that the large uncertainty 
in the mass estimate comes from the uncertainty in the tem- 
perature distribution and therefore part of the reason of the 
small scatter originates from the fact that the uncertainties are 
correlated. This is not true for the determination of the gas 
mass which does hardly depend on the temperature measure- 
ment. Therefore the Mgas.soo-T" relation, for which we found 
acoi - 0.16 for Ig(Mgas) for the REFLEX-DXL sample, is the 
most remarkably tight coiTelation. 

The scatter of the scaling relations of the REFLEX- 
DXL sample confirms the recent studies in observations (e.g. 
Reiprich & Bohringer 2002) and simulations (e.g. Borgani 
et al. 2004). For example, we obtained the coiTelative scat- 
ter of (0.19,0.20) for (lg(L),lg(M)) in the L-M500 relation and 
confirm the recent studies in observations (e.g. Reiprich & 
Bohringer 2002) and simulations (e.g. Borgani et al. 2004). 
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5. Discussion 

5. 1. Gas profiles in the outsl<irts 

The generally adopted yS-model (yS = 2/3) gives oc r~^. Mass 
in-falling becomes significant in the outskirts and thus makes 
the electron density profile steeper. Vikhlinin et al. (1999) 
found a mild trend for ji to increase as a function of cluster 
temperature, which gives ji ~ 0.80 and oc r"-^-^ for clus- 
ters around 10 keV. Bahcall (1999) also found that the electron 
number density scales as rig oc r'^-'^ at large radii. We confirm 
their conclusion that oc ^-2-42 at r > 2' for the REFLEX- 
DXL clusters. Due to the gradual change in the slope, one 
should be cautious to use a single slope double-yS model which 
might introduce a systematic error in the cluster mass measure- 
ments (as also described e.g. in Homer 2001). 

5.2. Validity of the spherical symmetry and hydrostatic 

equilibrium 

The total mass can be underestimated due to the assumption of 
spherical synmietry for elongated clusters (Castillo-Morales & 
Schindler 2003). However, a low surface brightness extension 
is often difficult to subtract correctly because of its less signif- 
icant boundary relative to the surroundings. In the REFLEX- 
DXL sample, RXCJ05 16.7-5430 is such an extreme case. A 
compression of the photon distribution extends from the clus- 
ter center to the north and a low surface brightness extends to 
rsoo south of the cluster center. We investigated this cluster and 
found the global measurements such as M500 are insensitive to 
the inclusion or removal of the region of the low surface bright- 
ness extension. The difference of the surface brightness profiles 
including and excluding the low surface brightness extension is 
within the l-cr error bar of the surface brightness. 

When X-ray images display a pronounced elongated 
and distorted cluster morphology, the cluster central po- 
sition cannot be determined unambiguously. For example, 
RXCJ0014.3-3022, RXCJ05 16.7-5430, RXCJ0528.9-3927, 
RXCJ0658.5-5556 and RXCJ1131.9-1955 are such cases in 
the REFLEX-DXL sample. For the case RXCJ0528.9-3927, 
we found that the measurements still agree with each other at 
the l-cr confidence level using the cluster center of the main 
component and the center of mass of the two components, re- 
spectively. 

Clusters showing significant merger features are still dy- 
namically young. These merger features can invaUdate the hy- 
pothesis of spherical symmetry and hydrostatic equilibrium. 
For example, RXCJ0658. 5-5556 is one of the most spectacular 
merger clusters. We checked the substructure excision method 
in this work by comparing the results including and exclud- 
ing the substructure in RXCJ0658.5-5556. We found that the 
global measurements are insensitive to the method excising 
substructure within the l-cr observational error. 

5.3. RXCJ0658.5-5556 and the scaling relations 

It is worthy to take a closer look, how clusters of 
particular morphological type affect the scaling rela- 



tions. We find that the pronounced merger clusters 
(RXCJ0014.3-3022, RXCJ05 16.7-5430, RXCJ2337.6-I-0016 
and RXCJ0658. 5-5556), which also show an excess in the 
central entropies and/or large gas mass fractions, introduce a 
significant broadening of the scatter in the scaling relations. 
Therefore the scatter has to be considered with caution since 
it could partially be an artificial effect of the invalidity of 
hydrostatic equiUbrium. A particular case is the extremely 
hot cluster, RXCJ0658. 5-5556. For example, excluding 
RXCJ0658. 5-5556 the Mgas,5oo-r relation provides reduced 
scatter, o'cor = 0.11 and ^cor = 0.06 for Mgas,5oo and T, 
respectively. 

RXCJ0658. 5-5556 has a very large weight on the slope of 
the scaUng relations due to its extreme location on the param- 
eter scale. It is also one of the merger systems with large un- 
certainties of the observational parameters. Therefore we have 
to be careful with the slope fitting when such merger clusters 
are involved. As an extreme case, the 5-7" relation of all the 
REFLEX-DXL clusters gives a slope twice higher than the 
pubUshed empirical slope. In this work, we either fitted the 
slope parameters excluding merger clusters (e.g. S-T) or fixed 
the slope parameters to those in the pubhshed papers (e.g. M- 
T). 

Merger clusters might also affect the determination 
of the normalization of the scaling relations. Excluding 
RXCJ0658.5-5556, the normaUzation of the scaUng relations 
of M50Q-T, Mgas,500-7", and l^T will be reduced by 5%, 5%, 
and 4%, respectively. RXCJ0658. 5-5556 lays above the nor- 
maUzed relations by a factor of 1.9,1.5, and 1.1 for Msoo-T", 
,500-7", and L-7', respectively. 

5.4. Dense core or cool gas 

A dense gaseous cluster core, as observed in the CCCs, does 
not necessarily require a cusp of the DM distribution as de- 
scribed by the extended-NFW model. Alternatively, a suffi- 
ciently cooler central temperature also results in a dense core 
without a cuspy DM profile which is demonstrated as follows. 
Restricting the analysis to r < 0.5', we use the pronounced 
CCC, RXCJ0307.0-2840, to illustrate the total mass distribu- 
tion in the dense and cool gaseous region in the cluster. In gen- 
eral, no cool gas has been observed showing a central tempera- 
ture lower than half of the mean temperature. We thus assume 
an extreme temperature drop to 1/3 of the observed temperature 
of the inner most bin towards the center. Assuming hydrostatic 
equilibrium, we derived a relatively low mass concentration in 
the cluster center with the steep central gas density distribution. 
We can not easily distinguish between the above two cases us- 
ing the current observations. 

5.5. Additional physical processes 

We observed a deviation around the self-similar model in the 
central region in the scaled profiles of the temperature, surface 
brightness, entropy, total mass and gas mass. This is most prob- 
ably the effect of different physical processes rather than sim- 
ply being statistical fluctuations in the measurements (Zhang 
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et al. 2004b, 2005b; Finoguenov et al. 2005). Many studies 
(e.g. Markevitch et al. 2002; Randall et al. 2002; Finoguenov et 
al. 2005) show that the X-ray property estimates in the center 
can be biased by phenomena such as ghost cavities, bubbles, 
shock and cold fronts, that may somehow invalidate the hydro- 
static equilibrium hypothesis and the assumption of homoge- 
neous temperature and density distributions. Complex dynami- 
cal interactions with AGN activities have been indicated by the 
coincidence of CCCs and radio sources (Clarke et al. 2005). 
M87 shows an example to test the effect of AGN interaction 
on the X-ray luminosity and also multi-temperature structure 
(Matsushita et al. 2002). 

After subtraction of the effect of measurement uncertain- 
ties, the remaining intrinsic scatter in the scaling relations is 
a signature in variations of cluster structure and ICM pro- 
cesses. This also includes merging clusters as an extreme case 
(Zhang et al. 2004b). Systematic studies of X-ray mergers have 
been done in observations using a series of cluster samples 
(e.g. Schuecker et al. 2001) and simulations (e.g. Schindler 
& Mueller 1993). Such a detailed study of the REFLEX- 
DXL clusters using a 2-dimensional approach can be found in 
Finoguenov et al. (2005). 

6. Summary and conclusions 

X-ray luminous (massive) clusters are used in a variety of 
ways to perform both cosmological and astrophysical studies. 
We selected an unbiased, almost volume-limited sample, the 
REFLEX-DXL cluster sample, from the REFLEX survey. We 
performed a systematic analysis to measure the X-ray observ- 
ables based on high quality XMM-Newton observations, and 
investigated various X-ray properties and the scaling relations 
of the REFLEX-DXL cluster sample. We summarize two main 
conclusions as follows. 

(i) An almost self-similar behaviour of the scaled profiles 
of X-ray properties, such as temperature, surface brightness, 
entropy, gravitational mass, and gas mass has been found above 
0.2-0.3 rsoo for the REFLEX-DXL sample. 

- The average global metallicity is 0.24 ± 0.032© for the 
whole sample. No significant evolution was found up to 
z ~ 0.3 in the metallicity comparing the REFLEX-DXL 
sample to the nearby galaxy cluster samples. This agrees 
with the results in Allen & Fabian (1998), Tozzi et al. 
(2003), Chen et al. (2003), De Grandi et al. (2004) and 
Pointecouteau et al. (2004). The results fit into the scenario 
showing no significant evolution of the iron abundance up 
toz ~ 1.1. 

- Based on the XMM-Newton observations, we obtained an 
average temperature profile of the REFLEX-DXL clusters, 
which agrees with the previous studies within the observa- 
tional dispersion. Markevitch et al. (1998) found a steep 
temperature drop beyond an isothermal center based on 
ASCA data. De Grandi & Molendi (2002) derived a uni- 
versal temperature profile which shows a similar decline 
using BeppoSAX observations. Based on high resolution 
Chandra observations, VikhUnin et al. (2005) and Piffaretti 



et al. (2005) confirmed the previous studies within the 
observational dispersion, where they found a more pro- 
nounced drop outside of 0.2-0.3 r2oo in the mean of the 
universal temperature profile. Additionally, Borgani et al. 
(2004) reproduced a similar temperature profile in their 
simulations as found by Vikhlinin et al. (2005). For the 
REFLEX-DXL sample, we found a universal temperature 
profile with a peak around 0.2-0.3 rsoo. We observed an 
almost constant value up to 0.3r5oo and a very mild de- 
crease to 80% of the peak value at r > 0.3r5oo for the 
average temperature radial profile of the REFLEX-DXL 
clusters. However, three clusters (RXCJ0437. 1+0043, 
RXCJ0658.5-H5556, and RXCJ2308.3+0211) show tem- 
peratures rising with relatively large error bars in their tem- 
perature profiles. The average temperature profile declines 
to 50% of the peak value at radii above 0.3r5oo when these 
three clusters are excluded. No significant detection of cool 
gas has been observed showing a central temperature lower 
than half of the mean temperature. 

- We determined the XMM-Newton surface brightness pro- 
files of the REFLEX-DXL clusters up to at least rsoo. We 
observe steeper profiles at large radii than what is gener- 
ally obtained for the j0-models with /3 ~ 2/3. The surface 
brightness profiles of most REFLEX-DXL clusters show a 
flat core populating a broad range of values up to 0.2 r^oo. 
The total luminosity of the cluster core (< 0.1 rsoo) ac- 
counts for about 8-33% cluster luminosity. For the clus- 
ters showing cooling cores, no well-defined constant cen- 
tral density was observed with the XMM-Newton resolu- 
tion. However, the surface brightness profiles are quite self- 
similar at r > 0.2r5oo for the REFLEX-DXL sample. 

- We performed the redshift evolution correction to the en- 
tropy profiles at 0.1 rsoo for the REFLEX-DXL clusters and 
obtained consistency with those for the nearby clusters in 
Ponman et al. (2003) (see Fig.|4]and Fig.fTOli within the ob- 
servational dispersion. However, four clusters showing rel- 
atively cool cores introduce a significant scatter in terms of 
a lower (~ 7%) normalization of the S-T relation, and three 
merger clusters introduce some scatter in terms of a higher 
(~ 10%) normalization in the S-T relation. Excluding 
these 7 clusters, the REFLEX-DXL sample shows an em- 
pirical scaling, 7'0 63±o.i5 shown in Fig. [10] the entropy 
profiles at r > 0.1 rsoo show a similar slope {S oc ^0-99±o.{)6-) 
as observed in Ettori et al. (2002b), S oc r" ''^. At the 1-cr 
confidence level, this scaling is shallower than the predic- 
tion of the spherical accretion shock model, S oc r' ' (e.g. 
Tozzi & Norman 2001; Kay 2004). We found that merger 
clusters show high central entropies, and the relatively re- 
laxed clusters show lower central entropies. Therefore the 
central entropy can be used as one means to distinguish the 
cluster dynamic state. For example, the observational de- 
viation of the (S ,T) data pair from the S-T scaling rela- 
tion for the individual cluster can be used to distinguish the 
relaxation stage. This may be particularly useful to char- 
acterize cluster merging along the line-of-sight, where the 
merger morphology is less obvious in the projection on the 
sky. 



12 



Zhang et al.: X-ray properties in massive galaxy clusters 



- We found a self-similar gravitational mass distribution for 
the REFLEX-DXL sample at radii above 0.2-0.3 rsoo. The 
precise gas density and temperature radial profiles provide 
a detailed diagnostics of the cluster structure and yield reli- 
able determinations of the total mass and gas mass fraction. 
In the outskirts, the observational density distribution pro- 
vides an average of p(r) oc 42 for the REFLEX-DXL 
sample. 

(ii) The scaling relations of the REFLEX-DXL sample at 
z ~ 0.3 agree with the scaling relations of the nearby samples 
within the observational dispersion after the redshift evolution 
correction. 

- Since the cluster temperatures of the REFLEX-DXL sam- 
ple are in a narrow temperature range, the whole sample 
provides a good means to constrain the normalization of 
the scaling relations (see Table |6]l and to study their intrin- 
sic scatter (see Table|7]l at the high mass end. 

- The results for the scaling relations of the REFLEX- 
DXL sample show good agreement compared with previ- 
ous studies (Table |6]l. This fits the general opinion (e.g. 
Maughan et al. 2003; Arnaud 2005; Vikhlinin et al. 2006) 
that the evolution of galaxy clusters up to z ~ 1 is well 
described by a self-similar model for massive clusters. 

- We found that the scatter of the normalization of the scal- 
ing relations is very sensitive to the cluster morphology. For 
example, the scatter of the L-r relation is reduced by 15% 
excluding 4 REFLEX-DXL clusters showing pronounced 
radiative cooling. Also the normalization of the L-T rela- 
tion excluding those 4 clusters is reduced by 10% for the 
REFLEX-DXL clusters. We investigated the logarithmic 
intrinsic scatter of the scaling relations which, for exam- 
ple, gives (0.19,0.20) for (L,M) in the L-Msoo relation and 
confirms the recent studies, such as (0.29,0.22) in Reiprich 
& Bohringer (2002). 
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Table 1. Properties of the suspicious point-like sources in the cluster center. Column (1): cluster name; Cols. (2,3): sky coordi- 
nates in epoch J2000 of the point-Uke source; Col. (4): index of the best power law fit; Col. (5): X-ray luminosity of the point-like 
source. 



RXCJ 


X-ray 


centroid 


Power law 


Lx(0. 1-2.4 keV) 




R.A. 


decl. 




lO'^ erg S-' 


0232.2-4420 


02 32 18.6 


-44 20 48.2 


1.77 ± 0.03 


5.22 ± 0.27 


0437.1-H0043 


04 37 09.8 


-hOO 43 48.9 


1.82 ±0.03 


3.56 ±0.18 


0528.9-3927 


05 28 52.6 


-39 28 16.8 


1.70 ±0.05 


3.20 ± 0.27 


1131.9-1955 


11 31 54.2 


-19 55 39.8 


1.75 ± 0.04 


1.97 ±0.08 


2308.34-0211 


23 08 21.6 


-02 11 29.1 


1.67 ±0.05 


2.24 ±0.12 


2337.6-1-0016 


23 37 35.3 


-1-00 15 52.1 


1.80 + 0.06 


1.02 + 0.10 



Table 2. Classification of the dynamical state. Column(l): classification; Col. (2): cluster names. 



Classification 


RXCJ 












Single 


0307.0- 


-2840, 


0532.9- 


-3701, 2308.3- 


-0211 




Primary with small secondary 


0232.2- 


-4420, 


0303.7- 


-7752 






Elliptical 


0043.4- 


-2037, 


0437.1+0043, 0516.7- 


-5430, 


1131.9-1955 


Offset center 


0014.3- 


-3022, 


0528.9- 


-3927, 0658.5- 


-5556, 


2337.6+0016 


Complex 


2011.3- 


-5725 











Classifications see Jones & Forman (1992). 

Table 3. Primary parameters. Column (1): cluster name; Col. (2): optical redshift (Bohringer et al. 2004); Cols. (3,4): sky coor- 
dinates in epoch J2000 of the cluster center; Col. (5): hydrogen column density in units of lO^'^cm"^ (Dickey & Lockman 1990); 
Col. (6):X-ray measured redshift; Col. (7): truncation radius corresponding to a of 3 of the observational surface bright- 
ness profile; Cols. (8,9): cluster global temperature and global metallicity in the 0.1-0.5 r^oQ region; Cols. (10,11): bolometric 
luminosity including and excluding cluster core (< O.lrsoo). 



RXCJ 


^opt 


X-ray centroid 


Nh 


^X-ray 


n 


T 


Z 


T incc 

Hoi 


rexcc 
Hoi 






R.A. 


decl. 


10^0 cm-2 




arcmin 


keV 


Zo 


10« 


erg s"' 


0014.3-3022 


0.3066 


00 14 18.6 


-30 23 15.4 


1.60 


0.28 ± 0.01 


7.61 


10.1 ± 0.3 


0.17 ±0.05 


2.12 ±0.17 


1.91 ± 0.16 


0043.4-2037 


0.2924 


00 43 24.5 


-20 37 31.2 


1.54 


0.28 ± 0.01 


7.02 


7.7 ± 0.3 


0.11 ±0.06 


1.07 ±0.11 


0.89 ± 0.09 


0232.2-4420 


0.2836 


02 32 18.8 


-44 20 51.9 


2.49 


0.27 ± 0.01 


6.81 


7.0 ± 0.3 


0.10 ±0.05 


1.89 ±0.14 


1.27 ±0.11 


0303.7-7752 


0.2742 


03 03 47.2 


-77 52 39.0 


8.73 


0.26 ± 0.01 


7.53 


8.2 ±0.5 


0.32 ± 0.08 


1.29 ±0.13 


1.04 ± 0.11 


0307.0-2840 


0.2578 


03 07 02.2 


-28 39 55.2 


1.36 


0.24 ± 0.01 


4.70 


6.4 ±0.3 


0.13 ± 0.06 


1.34 ± 0.13 


1.00 ± 0.10 


0437.1+0043 


0.2842 


04 37 09.5 


+00 43 54.5 


8.68 


0.31 ±0.02 


7.73 


5.1 ±0.3 


0.30 ±0.11 


0.62 ± 0.07 


0.42 ± 0.06 


0516.7-5430 


0.2943 


05 16 35.2 


-54 30 36.8 


6.86 


0.28 ± 0.01 


6.39 


7.5 ± 0.3 


0.16 ±0.07 


0.92 ±0.12 


0.84 ±0.11 


0528.9-3927 


0.2839 


05 28 52.5 


-39 28 16.7 


2.12 


0.26 ± 0.02 


5.80 


7.2 ± 0.4 


0.42 ± 0.08 


1.58 ±0.16 


1.25 ±0.13 


0532.9-3701 


0.2747 


05 32 55.9 


-37 01 34.5 


2.90 


0.27 ± 0.01 


7.01 


9.5 ± 0.4 


0.32 ± 0.07 


1.35 ± 0.12 


0.95 ±0.10 


0658.5-5556 


0.2965 


06 58 30.2 


-55 56 33.7 


6.53 


0.29 ± 0.01 


8.88 


10.6 ± 0.2 


0.26 ± 0.03 


4.87 ± 0.24 


4.21 ±0.21 


1131.9-1955 


0.3075 


11 31 54.7 


-19 55 40.5 


4.50 


0.33 ± 0.03 


9.00 


9.2 ± 0.4 


0.27 ± 0.05 


1.80 ±0.15 


1.58 ±0.14 


2308.3-0211 


0.2966 


23 08 22.3 


-02 11 32.1 


4.45 


0.27 ± 0.02 


7.31 


7.9 ± 0.7 


0.31 ±0.14 


1.20 + 0.13 


0.88 + 0.11 


2337.6+0016 


0.2753 


23 37 37.8 


+00 16 15.5 


3.82 


0.27 ± 0.01 


7.27 


9.6 ± 0.3 


0.25 ± 0.06 


1.00 ± 0.09 


0.79 ± 0.08 
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Table 4. Deduced properties of the REFLEX-DXL clusters. Column (1): cluster name; Col. (2): central electron number density; 
Col. (3): central entropy; Cols. (4,5): cooling time and cooling radius; Col. (6): rsoo; Cols. (7-9): gas mass, total mass, and gas 
mass fraction at r^oo- 



RXCJ 




J 




'cool 


'' 500 


■'*^ga.s,500 




7gas,500 




CI 0-3 rm-^^ 


CkeV cm^ 














0014 'i-T.no'y 


<i 7 + n 9 


491 + <17 


14.1 + 0.6 


00 


1.24 


1.0 + 0.2 


/ .*+ in Z..7 


1 49 + 094 


004^ 4—90^7 


5.8 + 0.4 


293 + 17 




0.12 


1.08 


0.6 + 0.1 


4 R + 1 R 


190 + 09^ 


09^9 l—AAOO 

\JjLJjL,jL — ' 1 ' 1 jlAJ 


14.0 + 0.1 


970 + in 


10 3 + 09 




1 30 


Q + n 9 


8.4 + 2.4 


lOS + 09^ 


VJJVJJ, III JZ, 




^11 +94 


10 3 + 04 


19 


1 97 


8-1-01 

U.O ^ U. 1 


7 7 j_ 9 'J 


101 -1-0 099 


0307.0-2840 


10.2 ± 0.3 


265 ± 16 


10.5 ± 0.3 


0.17 


1.14 


0.6 ±0.1 


5.5 ± 1.1 


0.113 ± 0.029 


0437.1+0043 


10.1 ±0.4 


223 ± 14 


10.2 ± 0.3 


0.16 


1.17 


0.5 ±0.1 


6.1 ± 2.2 


0.081 ± 0.020 


0516.7-5430 


2.9 ± 0.2 


403 ± 24 


16.7 ± 0.4 


0.00 


1.19 


0.8 ± 0.2 


6.4 ± 2.1 


0.122 ±0.029 


0528.9-3927 


10.6 ± 0.4 


280 ± 16 


10.3 ± 0.3 


0.14 


1.19 


0.9 ±0.1 


6.4 ± 1.5 


0.135 ±0.028 


0532.9-3701 


13.3 ±0.5 


337 ± 23 


10.4 ±0.3 


0.14 


1.13 


0.6 ±0.1 


5.4 ± 1.7 


0.107 ±0.024 


0658.5-5556 


6.1 ±0.2 


354 ± 16 


10.1 ±0.2 


0.15 


1.42 


1.8 ±0.3 


11.0±6.8 


0.161 ±0.018 


1131.9-1955 


5.8 ±0.3 


307 ± 18 


10.0 ± 0.3 


0.11 


1.10 


0.8 ± 0.2 


5.2 ±3.0 


0.160 ±0.030 


2308.3-0211 


9.2 ± 0.3 


296 ± 22 


10.1 ± 0.3 


0.14 


1.24 


0.7 ±0.1 


7.4 ± 1.8 


0.089 ± 0.022 


2337.6+0016 


5.7 ±0.3 


372 ± 21 


12.3 ±0.3 


0.00 


1.43 


0.8 ±0.1 


10.9 ± 2.6 


0.074 ± 0.016 


Mean 
















0.116 ±0.007 



Table 5. Parameters of the extended-NFW model fits of the REFLEX-DXL clusters. Column (1): cluster name; Cols. (2,3): 
characteristic density and scale of the halo; Cols. (4): slope parameter; Col. (5): and its degree of freedom (d.o.f.); Col. (6): 
concentration parameter at ^500. 



RXCJ 


Ps 

(Me 


Mpc-3) 




''s 

(Mpc) 


or 


ArVd.o.f. 






C500 


0014.3-3022 


(2.7 


± 0.3) X 


10^ 


0.16 ±0.01 


-1.69 


90.0/157 


7, 


.5 


± 1.5 


0043.4-2037 


(1.8 


± 0.2) X 


lO"* 


0.15 ± 0.01 


-0.56 


30.9/115 


7, 


.4 


± 1.5 


0232.2-4420 


(1.3 


± 0.7) X 


10^ 


1.13 ± 0.74 


1.28 


65.7/124 


1, 


.2 


± 0.4 


0303.7-7752 


(8.8 


± 0.6) X 


10' 


0.22 ± 0.01 


-0.17 


18.0/143 


5, 


.7 


± 1.2 


0307.0-2840 


(3.7 


± 0.2) X 


10" 


0.11 ±0.01 


-0.87 


3.8/82 


10, 


.7 


±2.2 


0437.1+0043 


(6.2 


± 0.7) X 


103 


0.23 ± 0.01 


0.08 


34.8/97 


5, 


.1 


± 1.1 


0516.7-5430 


(2.5 


± 0.2) X 


103 


0.38 ± 0.02 


-0.09 


21.4/131 


3, 


.1 


±0.6 


0528.9-3927 


(4.6 


± 2.6) X 


10"* 


0.08 ± 0.02 


-0.51 


113.8/56 


14, 


.6 


±4.4 


0532.9-3701 


(6.9 


± 1.1) X 


10" 


0.80 ± 0.01 


-0.65 


39.7/118 


14, 


.2 


±3.0 


0658.5-5556 


(4.0 


± 0.6) X 


103 


0.43 ± 0.03 


-0.27 


333.6/240 


3 


.3 


±0.7 


1131.9-1955 


(1.2 


± 0.3) X 


10« 


0.04 ± 0.01 


-7.71 


49.9/128 


30, 


.0 


±6.6 


2011.3-5725 


(4.3 


±4.1)x 


10* 


0.07 ± 0.03 


-0.44 


7.0/8 


4, 


.2 


±0.5 


2308.3-0211 


(2.4 


± 0.3) X 


10* 


0.13 ± 0.01 


-0.31 


9.4/58 


9 


.7 


± 2.0 


2337.6+0016 


(6.8 


± 1 .9) X 


10' 


0.23 ± 0.03 


0.12 


82.1/91 


6 




± 1.5 
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X 



Reference 



^ E{z) (A,,JA,.o)°-^ 


T 
keV 


lO^^s/jji 


1.5 (fixed) 


Evrard et al. 96 






io"-«-§:g?/,jgi 


1 

i-^o_0.12 


Finoguenov et al. 01b 






10l3.96±0.02;j-^l 


1.5 (fixed) 


Ettori et al. 04 






1013.46±0.02;,-^1 


1.59 ± 0.05 


Popesso et al. 05 






ir|13.57+0.02;,-l 


1.49 ±0.15 


Amaud et al. 05 








1.5 (fixed) 


This work 



"185,500 r\ ;a ^0.5 
Mo (A^/AcO) 


T 
keV 




1fll2.63±0.03l,-l 


1 QQ _i_ n 1 Q 
l.Vo + U.io 


Monr et al. yy 






lU "50 
inl2.53±0.04l,~l 

iU /2yQ 


1.80 ±0.16 

1 c ^'f;vorl^ 
i.o (^nxeaj 


Castillo-Morales & Schindler 03 
This work 


rU.I-i.4kcV 
ergs-1 (Ac,z/A,,o) "-^ 


r 




1 r\49 04 » -1 

lU «JQQ 


2.10 ± 0.24 


Markevitch 98 






iU ftjo 
1042.19±0.29;,-^1 

1042.37±0.06;,-^1 


Z.OU ± U. i J 

2.44 ± 0.39 
2.60 (fixed) 


Reiprich & Bohringer 02 
Ikebe et al. 02 
This work 


^ £(z)-i (A„/A,,o)-°-= 


T 
keV 




ir.4243±0.04;^-l 
1042.82±0.03;j-^l 
1042.97±0.07;j-^l 
1042.85±0.09;j-^l 
1042.38±0.06;,-^1 


2.64 ± 0.27 
2.88 + 0.15 
2.79 ± 0.08 
2.98 ±0.12 
2.98 (fixed) 


Markevitch 98 
Amaud & Evrard 99 
Xue & Wu 00 
Reiprich & Bohringer 02 
This work 


£(z)-' (A,,,/A,,o)-''-5 


Mo 


E{z) (AcJAcof 


1q22455±1.298;,-^1 


1.504 ±0.089 


Reiprich & Bohringer 02 




1025.19±0.10;,-^1 
1025.00±0.06;,-^1 


1.30 ±0.12 
1.3 (fixed) 


Popesso et al. 05 
This work 


^ E(z)-' (A„/A,o)-o-5 


Msoo 


Eiz) (A.JA.o)"' 


1025.35±0.06;,-^1 


1.3 (fixed) 


This work 








0.035 ± ().()30/)^J - 


0.5 ± 0.3 


This work 



Table 7. Logarithmic intrinsic scatter measured for the scaUng relations. Column(l): scaUng relation; Col. (2): variable; 
Cols. (3,4): logarithmic intrinsic scatter and correlative scatter. 



Rclalion 


Variable 


""inl 




Msoo-r 


lg(M500) 


0.17 


0.20 




m 


0.09 


0.10 


^gas,500-7' 


lg(^gas,5O0) 


0.14 


0.16 




m 


0.07 


0.08 


2^0.1-2.4keV_j' 


lg(^0.1-2.4keV) 


0.16 


0.17 




ig(r) 


0.06 


0.07 


r0.1-2.4keV 14 


lg(^ai-24keV) 


0.19 


0.20 




lg(M5oo) 


0.16 


0.19 




lg(ibol 


0.16 


0.17 




lg(T) 


0.06 


0.07 


Lbol-Msoo 


lg(ibol 


0.18 


0.19 




lg(M5oo) 


0.17 


0.20 
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_j . . . . i ■ ■. , I •rn. a . . 

0.5 1 5 10 

channel energy (keV) 



Fig,l. XMM-Newton spectra (pn in grey and MOS in black) of RXCJ0014.3-3022 of the 0.1-0.5 r^oo region fitted by the 
combined "wabs*mekal" model. The ratios of the observational data to the models are in the lower parts of the panels (offset 
zero for pn, +1 for MOSl, +2 for MOS2). 
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Fig. 2. De-projected temperature profiles of the REFLEX-DXL clusters. The temperature profiles are approximated by the pa- 
rameterization T{r) = T-i exp[-(r - Tif-jT^] + T(,{\ + r^/T]^)''''^ + T-j crossing all the data points (solid). The l-cr confidence 
intervals are approximated by T{r)\\ ± {T% + Tgr)] (dashed). The vertical Une denotes O.lrsoo. 
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R (arcmin) R (aromin) 



Fig. 3. Surface brightness profiles of the REFLEX-DXL clusters and their best fits. Residuals scaled by the data uncertainties are 
plotted in the upper part of the diagrams. 
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RXCJ0014.3 
RXCJ0043.4 
RXCJ0232,2 
RXCJ0303.7 
RXCJ0307.0 
RXCJ0437.1 



3022 
2037 
4420 
7752 
2840 
+0043 



RXCJ0528.9-3927 
RXCJ0532. 9-3701 
RXCJ0658.5-5556 
RXCJ1 131.9-1955 
RXCJ2308.3-021 1 
RXCJ2337.6+0016 



T (keV) 



Fig. 4. Redshift corrected entropy (diamonds) at O.lrsoo vs. T for the REFLEX-DXL clusters. Four clusters showing relatively 
cool cluster cores are in thick lines. Nearby clusters of Ponman et al. (2003, boxes) are shown for comparison. The solid line 
denotes the best fits {S oc fO-6i±o.i5y electronic edition of the Journal for a colour version. 
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0.1 1 ~ 0.1 1 



Fig. 5. Measured mass profiles of the REFLEX-DXL clusters and their best fits by the extended-NFW models. Residuals scaled 
by the data uncertainties are plotted in the upper part of the diagrams. 
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Fig. 6. The concentration parameter of the extended-NFW model versus mass for the REFLEX-DXL sample. The line denotes 
the best fit of the REFLEX-DXL sample. The colours have the same meaning as those in Fig.H] 



2000 1000 500 200 

A 



Fig. 7. Gas mass fraction as a function of density contrast for the REFLEX-DXL sample. The l-tr error corridor (dashed curves) 
is displayed for RXCJ0532. 9-3701 only. The colours have the same meaning as those in Fig.H] 



This work, hatched 
Vikhlinir et ol. 05. filled 
Morkevitch et al. 98, hatched 




Fig. 8. Left: Scaled radial temperature profiles. The colours have the same meaning as those in Fig. |4l Right: An average temper- 
ature profile of the REFLEX-DXL clusters compared to the temperature profile ranges in Markevitch et al. (1998, grey, hatched) 
and Vikhlinin et al. (2005, grey, filled). 
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0.01 0.1 



Fig. 9. Surface brightness profile fits scaled according to the empirical scaling, Sx ^'^^ given in Arnaud et al. (2002). The 
colours have the same meaning as those in Fig.H] 




Fig. 10. Left: Scaled entropy profiles and the best fit of the r > O.lrsoo region. The merger clusters (RXCJ0014. 3-3022, 
RXCJ05 16.7-5430 and RXCJ2337. 6+0016) are shown with asterisk symbols. Tio denotes temperatures in units of 10 keV. 
Right: Scaled entropy profile fits of the REFLEX-DXL sample compared to the Birmingham-CfA clusters (grey shadow) in a 
temperature range of 6-20 keV. The colours have the same meaning as those in Fig.|4] 




Fig. 11. Scaled total mass profiles (left) and gas mass profiles (right) for the REFLEX-DXL sample, with the 1-cr error corridor 
(dashed) displayed for RXCJ0532. 9-3701. The colours have the same meaning as those in Fig.|4] 
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L^KObAi (0.1-2.4 keV) (erg s"') 

Fig. 12. A comparison of the XMM-Newton luminosity and ROSAT luminosity taken from Bohringer et al. (2004). Filled and 
open circles correspond to the substructure and point-like source subtracted and un-subtracted ROSAT luminosity, respectively. 
The point-like sources and substructures are detected using the XMM-Newton observations. The colours have the same meaning 
as those in Fig.|4] 




T (keV) 

Fig. 13. Mass vs. temperature for the REFLEX-DXL sample and the best fit power law (black). The grey lines denote the best 
fits in Evrard et al. (1996, sohd), Ettori et al. (2004, dashed), Popesso et al. (2005, dash-dotted) and Arnaud et al. (2005, dotted). 
The colours have the same meaning as those in Fig. |4l 
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Mohr et al. 99 

Castillo— Morales Sl Schindler 03 
This work 

1 /' 









5 10 20 

T (keV) 



Fig. 14. Gas mass vs. cluster temperature and the best fit power law (black). The grey lines denote the best fits in Mohr et al. 
(1999, dashed) and Castillo-Morales & Schindler (2003, solid). The colours have the same meaning as those in Fig. HI 
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T (keV) 



Fig. 15. X-ray luminosity in the 0. 1-2.4 keV band vs. cluster temperature and the best fit power law (black). The grey lines denote 
the best fits in Markevitch (1998, solid), Ikebe et al. (2002, dashed), and Reiprich & Bohringer (2002, dash-dotted). The colours 
have the same meaning as those in Fig.|4] 
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T (keV) 



Fig. 16. X-ray bolometric luminosity vs. cluster temperature and the best fit power law (black). The grey lines denote the best fits 
in Markevitch (1998, solid), Arnaud & Evrard (1999, dashed), Xue & Wu (2000, dash-dotted), and Reiprich & Bohringer (2002, 
dotted). The colours have the same meaning as those in Fig.|4] 




Fig. 17. X-ray luminosity in the 0.1-2.4 keV band vs. total mass and the best fit power law (black). The grey lines denote the 
best fits in Reiprich & Bohringer (2002; dashed) and Popesso et al.(2005, solid). The colours have the same meaning as those in 
Fig.a 
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Fig. 18. X-ray bolometric luminosity vs. cluster total mass and the best fit power law (black). The colours have the same meaning 
as those in Fig.|4] 




T (keV) 



T (keV) 



Fig. 19. Gas mass fraction /gas,500 (left) and Gas mass fraction ratio /gas,5oo//gas,2500 (right) as a function of cluster temperature. 
The best fit is /gas,500 t'^-^'^'^-^. The colours have the same meaning as those in Fig.|4] 
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-0.05 0.05 

(^int(^^gos,50o) 



0.1 



Fig. 20. Histogram of the intrinsic scatter for lg(Mgas,5oo) of the Mgas.soo-T" relation for the REFLEX-DXL clusters. 



